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Abstract 

In this paper we establish a relation between the spread of infectious diseases and the dy- 
namics of so called M/G/1 queues with processor sharing. The in epidemiology well known 
relation between the spread of epidemics and branching processes and the in queueing theory 
well known relation between AZ/G/l queues and birth death processes will be combined to pro- 
vide a framework in which results from queueing theory can be used in epidemiology and vice 
versa. 

In particular, we consider the number of infectious individuals in a standard SIR epidemic 
model at the moment of the first detection of the epidemic, where infectious individuals are 
detected at a constant per capita rate. We use a result from the literature on queueing processes 
to show that this number of infectious individuals is geometrically distributed. 

1 Introduction 

1.1 Epidemiological motivation for the paper 

In real-life, the knowledge that an infectious disease is spreading may change the further spread 
of the disease. Both deliberate intervention measures and non-mandatory behavioural changes 
can contribute to this phenomenon [1]. Examples are the 2002-2003 SARS epidemic, during which 
people avoided crowded places; the nosocomial pathogen methicillin-resistant Staphylococcus aureus 
(MRSA), where in the Nordic countries and The Netherlands hospitalised patients known to be 
colonised with MRSA are treated in single-bed isolation rooms and all patients who might have 
had contact with the index case, i.e., at least all patients in the same hospital ward as the detected 
patient, are screened for colonization [2j; or contagious animal diseases like Foot and Mouth Disease 
or Classical Swine Fever, where farms are depopulated at the moment the disease is detected 1^. 
The distribution of the number of infectious individuals at the moment of the first detection of 
the epidemic is important. Not only because the disease dynamics before the first detection is not 
influenced by control measures and, therefore, can be a measure for the true spreading capacity 
of the disease but also because this distribution informs us on the amount of control required to 
curtail the outbreak. 

In this paper we deal with the change of the epidemic process at detection, by exploring an 
obvious but not so well-known relationship between epidemics and queueing theory. By using this 
relationship and results from queueing theory we are able to derive the distribution of the number 
of infectious individuals at the moment of the first detection in a broad class of epidemics in large 
populations. Furthermore, the approach of tackling problems from epidemiology by using queueing 
theory is promising itself and might lead to results on epidemics that are beyond the scope of this 
paper. 

1.2 The stochastic processes 

We first consider an SIR (Susceptible — > Infectious Removed/Recovered) epidemic [H |5] with 
detections in a fixed population of size n. In this model each pair of individuals contacts each other 
at a constant rate of A/(n — 1). This implies that every individual makes contacts at a total rate A, 
i.e., the number of pairs it belongs to times the contact rate per pair. If a contact is made between 
an infectious and a susceptible individual, then the susceptible individual becomes infectious. Note 
that we define contacts as events at which infection will take place if one of the individuals is 
infectious and the other is susceptible. These contacts need not be the same as physical contacts. 
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If the physical contact rate is A'/n and the probabihty of transmission during a physical contact is 
c, then A = cA'. 

After an infectious period which is distributed as the random variable L, and independent of 
other infectious periods, an infectious individual becomes removed, which means that the individual 
becomes immune and stays so forever. During its infectious period an individual might be detected, 
which happens at a per capita rate of 5. At the moment of a detection, the process might change. 
However, because we are interested in the number of infectious individuals at the moment of a 
detection, it is in the context of this paper not important what happens after the detection. 

This process can be related to a branching (birth-death) process [6], by interpreting an individual 
(say j) that is infected by another individual (say i) as a child of i. This interpretation relates the 
number of infectious individuals at time t to the number of individuals in the branching process at 
time t. For finite population size n, this is not a proper branching process, because the distribution 
of children per individual changes over time because of depletion of susceptibles. However, for 
large n, the probability that during the initial phase of the epidemic contacts are made between 
individuals that are both non-susceptible is small and, therefore, the start of an epidemic can be 
analysed using the corresponding branching process [HlZj. In this branching process individuals 
have independent random life lengths, all distributed as L. During their lifetime they get children 
at a constant rate A. As in the epidemic process, detections take place at a per capita rate of 5. 

Remark: For definition purposes we only consider the SIR epidemic in this paper. However, all 
results of this paper hold for all types of epidemics in which the number of infectious individuals, 
can be approximated by the number of living individuals in a branching process in which both 
the per capita birth rate and the per capita detection rate are constant over the lifetime. So, 
also the SIS (Susceptible — > Infectious — > Susceptible) epidemics, in which recovered individuals 
are immediately susceptible again, SIRS (Susceptible — > Infectious — > Removed/Recovered — > 
Susceptible) epidemics, in which recovered individuals become susceptible again after some (possibly 
random) time period or epidemics in homogeneous populations in which deaths and births are taken 
into account may be considered. 

From queueing theory, we know that we can relate the so-called M(emoryless)/G(eneral)/l- 
queue with a Processor Sharing service discipline (or M/G/l-PS) queue O [9l [lOl [11] to a branching 
process. The M/G/1 queue is defined as follows. Customers enter a queue at a constant rate A 
and they require a random amount of serving time (their workload) from a server. The workloads 
of the customers are independent and distributed as L. A single server is serving the queue. 

The processor sharing service discipline may be interpreted as the limit of a round robin service 
discipline. In that discipline, the server serves a customer either for a time length h, or if the 
remaining amount of serving time is less than h, until the customer is completely served. After 
this time length the server moves on to the next customer. Customers that are not completely 
served yet stay in the queue, while completely served customers leave the queue. After serving the 
last customer in the queue, the server returns to the first customer that is still in the queue. The 
processor sharing discipline is obtained by taking the limit /i \ 0. 

For the M/G/1 queue with round robin service discipline, we may consider customers that 
arrive during the time that the server was serving customer i as children of i. In this way we obtain 
a branching process, in which individuals can give birth to other individuals during a random 
period distributed as L and during which it gets children at a constant rate A. Note that time in 
the queueing process is not the same as time in the branching process, because in the branching 
process the individuals all grow older at a constant speed, while in the queueing process, customers 
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are served one by one, where the server has a constant speed. So, if the number of customers in 
the row is increasing, the service per time unit for a given individual in the queue is decreasing. 
However, if /i \ 0, the order of events (arrivals/births and departures/deaths) in the processes 
is the same. Another difference between the queueing process and the corresponding branching 
process is that for the branching process the state without individuals is absorbing, while in the 
queueing process customers may arrive in an empty queue, but up to the first time the queue 
becomes empty, the two processes are in direct correspondence. 

We add a "catastrophe" process [12] to the queueing process in which catastrophes occur at 
a constant rate 6. At the time of a catastrophe, a sudden change might happen to the queue, 
e.g., the queue might be emptied or a fraction of the queue might leave the queue. However, in 
this paper we are only interested in the number of customers in the queue at the moment of a 
catastrophe. Therefore, we do not need to specify what happens at the time of a catastrophe. Yet 
in the epidemic context we shall assume that detection has no impact at all, as this helps to keep 
the notation simple. 

The catastrophes are incorporated in the branching process corresponding to the M/G/ 1-PS 
queue as follows: Consider the M/G/1 queue with round robin service discipline. If a catastrophe 
occurs at a moment customer i is served, then it can be seen as a detection of individual i in the 
branching process, which occur at a per-capita rate 5. The M/G/l-PS queue with catastrophes 
is obtained by /i \ 0, where h is the time a server stays with the same customer. The order of 
events (arrivals/births, departures/deaths and catastrophes/detections) in the branching process 
with detections and the M/G/ 1-PS queue is still maintained. So, the number of individuals alive 
at the moment of the first detection in a branching process has the same distribution as the number 
of customers in the M/G/l-PS queue at the first catastrophe. 

We use a result from Kitaev [11], to show that the number of customers in a queue after an 
exponentially distributed time, conditioned on the queue being in its first busy period, i.e., the 
period in which the server was non-stop working, starting at the arrival of the first customer, is 
geometrically distributed. This implies that the number of infectious individuals at the moment of 
the first detection for an SIR epidemic in a large randomly mixing population is also geometrically 
distributed. So, the distribution of this number can be described by one parameter, while the 
process itself is described by 2 parameters and an unspecified distribution. This implies that if 
the only observations available, are the number of infectious individuals at the moment of the first 
detection in different outbreaks of the same disease, then we can only estimate 1 parameter and we 
cannot provide estimates for A, 6 and L separately. 

1.3 Earlier work 

The relation between M/G/1 queues and birth and death processes has already been discussed by 
Kendall in [QjITO]. In [11] a discussion on the relation between birth and death processes and the 
M/G/1 queues with processor sharing can be found. References to earlier work on this subject are 
also given in that paper. 

We are aware of only a few references to queueing theory in the epidemiological literature. In 
[7] results for the M/M/1 queue (the queue with exponentially distributed workloads) are used to 
compute the total cost of an epidemic, which is interpreted as a constant times the total number 
of infection days (i.e., it is interpreted as a constant times the integral of the number of infecteds 
over time). In [13] (see also [H P-12]) a construction for an SIR epidemic is given, which is very 
close to interpreting the epidemic as an M/G/1 queue with processor sharing, however the relation 
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is not made explicit there. The random time change discussed here was also applied in |14j . but in 
that paper was no explicit reference to queueing theory. 

1.4 Outline of the paper 

In the next section, we define the SIR epidemic with detections, the corresponding branching 
process with detections and the M/G/l-PS queue with catastrophes in mathematical language, in 
order to show the relationship between epidemics and queues in a rigorous way. In Section [3] we 
will provide a random time change argument to establish the relation between SIR epidemics and 
M/G/l-PS queues and we show that the number of infectious individuals at the moment of first 
detection in a large population is geometrically distributed. We use explicit computations or refer 
to literature in which the results are already proven rigorously. However, if possible, we also give 
heuristic and intuitive arguments for the claims made, which might be helpful for further use of 
queueing theory in epidemiology. In particular, in Section U] we provide an intuitive proof for the 
geometric distribution of the number of infectious individuals at the moment of first detection if 
the infectious periods are exponentially distributed. Some applications of the results of this paper 
are discussed in Section [5j In the final section we discuss possible extensions and some limitations 
of the use of queueing theory in epidemics. In particular, we discuss whether the coupling between 
branching processes and M/G/1 queues can be made for branching processes corresponding to 
epidemic models in which the infectivity/contact rate of an individual is not constant during its 
infectious period. 

2 Definitions and notation 

Before we formally describe the relevant processes for this paper, we give some definitions. We 
define 1/0 := lim 1/x = oo, 1/cxd := lim 1/x = and x cxd = lim x j; = 0. Furthermore, 

X'\0 X—fOO z— >oo 

the maximum/supremum of an empty set of real numbers is 0, while the minimum/infimum of 
an empty set is oo. The indicator function IL(A) takes value 1 if the event A occurs and if the 
event A does not occur. For a function f{x), we define f{x—) := lim /(y) and /(x) = o{x) if 

lim f{x)/x = 0. The natural numbers, not including 0, are denoted by N and Nq := N U {0}. 

x\0 

Finally, with some abuse of notation, for all processes under consideration, we will use {J-t] i > 0} 
to denote the filtration to which the process is adapted [15, p. 475]. So, all information on the 
process available at time t is contained in Tt- 

Throughout this paper we may deviate from standard notation in literature, because we want to 
relate different processes, which all have different standard notations for quantities that are related 
in this paper. In our notation we want to make clear which quantities in the different processes are 
related. 

2.1 The SIR epidemic with detections 

We consider an SIR epidemic in a homogeneous and randomly mixing population without demo- 
graphic turnover of size n. Let S^'^\t) be the number of susceptible individuals in the population 
at time t, I^'^\t) the number of infectious individuals at time t and R^'^\t) the number of re- 
moved/recovered individuals at this time. The epidemic starts with one infectious individual in a 
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further susceptible population, i.e., 5^"-'(0) = n — 1, /("^(O) = 1 and R^'^\0) = 0. Furthermore, we 
assume that the initial infectious individual was infected itself at time 0. 

Every pair of individuals makes contact at a constant, strictly positive, rate A/(n — 1), i.e., 
contacts between a pair of individuals are made according to a Poisson process with parameter 
A/(n — 1). If a contact is made between an infectious and a susceptible individual, the susceptible 
individual becomes immediately infectious. 

An infected individual stays infectious for a random infectious period, which is distributed as 
the random variable L, which is almost surely (a.s.) positive, i.e., P(0 < L < oo) = 1. The 
infectious periods are independent and identically distributed (i.i.d.). After the infectious period 
individuals become immune and stay so forever. 

We extend the standard SIR epidemic model by adding a detection process. In this process, 
infectious individuals are detected at a constant per capita rate 6. Mathematically, the detection 
process corresponds to a Poisson process with parameter 5 > 0, which is defined on all infectious 

(n) 

periods. Let denote the random time of the i-th detection in the population. If the number 

(n) 

of infectious individuals is before the i-th. detection takes place, we say = oo. As mentioned 
before, in this paper we are interested in the number of infectious individuals at the time of first 
detection I^^\D^^), conditioned on < oo. 

2.2 The branching process 

Let Z{t) := -^A,/3(s),<5(^) be a branching process p] with detections in which individuals have i.i.d. 
life lengths, distributed as L, with moment generating function /3(s) := E(e~'*^) and give birth 
at a constant per capita rate A > 0. So new individuals are born at a total rate of XZ{t). The 
process starts with one individual, whose age at time is 0. There is a detection process on top of 
this process of births and deaths, which is a homogeneous Poisson process with rate 5 > 0, which 
is defined on the life times of the individuals. So, detections happen at a total non-constant rate 
6Z{t). Again, detections do not influence the further course of the branching process. We only use 
the subscripts in ^a,/3(s),(5(^) if we want to stress the dependence on these parameters. 

We use Di to denote the time of the i-th detection in real-time in the branching process. If the 
branching process goes extinct before the i-th detection, then Di = oo. We label the individuals in 
the epidemic by the real-time order in which they appear in the branching process, so the ancestor 
in the branching process gets label 1, its first child gets label 2, then the next individual to be 
born (which is either a child or grand-child of the ancestor) gets label 3, and so on. Note that in 
general, it is not individual i that is detected at time D^. The random life length of individual i 
is denoted by Li and its time of birth by Tj. So, Ti = 0. We define Ai{t), the age of individual 



i at time t, as Ai{t) = {t — Ti)l(Tj < t) for t < Tj + Li and Ai{t) = oo if t > Tj + Lj. So, 



i=l 

is defined in terms of events occurring no later than t, so Z{t) G J-'t- 

For i £ N, the random variables Ti and Di are such that both of the sequences are increasing 



Ai{t) =1(0 < Aiit) < Li). 



Observe that 



oo 




(1) 
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in i and 

oo 

¥{Tie {t,t + h)\J^t,Z{t) = k,^l{Tj <t) = i-l) = Xkh + o{h), 

i=i 

oo 

¥{Di e {t,t + h)\J^t,Z{t) = k,^l{Dj <t) = 1-1) = Skh + o{h), (2) 

oo oo 

¥{^i{Tie{t,t + h)) + Y,'^{Die{t,t + h))>i) = o{h). 

i=l 1=1 
oo oo 

If Z{t) = and ^^Tj <t)=i (resp. ^%Dj <t) = i), then Tjt = oo (resp. Dk = oo) for k > i. 

j=i j=i 

2.3 The M{X)/G{l3{s))/l queue with Processor Sharing and catastrophes 

We define the M(A)/G(/3(s))/l queue with Processor Sharing with catastrophes or M/G/l-PS 
queue with catastrophes, Q{t) := Q (t) := ^^^^(s) 5^^)' ^ follows. Customers arrive in a queue 
according to a homogeneous Poisson process on the positive half line (0, oo), with rate A. Each 
customer brings in a workload, which is independent of the arrival process and workloads of other 
customers and distributed as the random variable L, with moment generating function /3(s). We 
assume that P(0 < L < oo) = 1. The customers in the queue are labelled according to the order in 
which they arrived in the queue. The time of arrival of the i-ih customer is denoted by T,. Unless 
specified otherwise, we assume that no customers are in the queue at time t = 0. 

The total workload customer i brings in is denoted by Li. The amount of service time customer 
i already received at time t, is denoted by Ai{t), where we define Ai{t) = for t < Tj. Customer 
i leaves the queue at sup{t;Ai{t) < Li} and from that time on (including the time of departure) 
Ai{t) = oo. 

One server serves the people in the queue in such a way that all customers in the queue receive 
the same amount of service per time unit. Thus, ^Ai(t) =1(0 < Ai{t) < oo)/Q{t), where Q{t) is 
the number of customers in the queue at time t. 

Independently of the ordinary M/G/l-PS queueing process as defined above, we define a catas- 
trophe process. Catastrophes occur according to a Poisson process on (0, oo) with rate 6. The time 
of the i-th catastrophe is denoted by Di. 

Note that 

oo 

Q{t) = J^mfi < t) -liMt) = oo)). (3) 

i=l 

So, Q{t) is (as is Z{t)) defined in terms of events occurring no later than t, i.e., Q{t) G J^f For 
i G N, the random variables Tj and Di are defined such that both of the sequences are increasing 
in i and 

oo 

F{fie{t,t + h)\J^t,^ii{fj <t) = i-i) = \h + o{h), 

oo 

¥{Die it,t + h)\J^t,y^ll{Dj <t) =1-1) = Sh + o{h), (4) 



oo 



¥{Y^i{fie{t,t + h)) + J2KDieit,t + h))>l) = o{h). 



i=l i=l 
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Note that, while Z(0) / 0, Q(0) = 0. 



3 The relationship between the SIR-epidemic with detections and 
the M/G/l-PS queue with catastrophes 

In this section we show that for the SIR epidemic with infection rate A, detection rate 6 and 
moment generating function of the infectious period (3{s) = E(e~'^^), 

lim P(/(")pi")) = klDi"^ < oo) = P(Z;,,^(,),5(Z)i) = fclTi = 0,Di< oo). 

n— >oo ^ ' 

So, in the large population limit, the number of infected individuals at the time of first detection 
is distributed as the number of alive individuals at the time of first detection in the corresponding 
branching process. 

After that we prove that 

nzx,mADi) = m = 0, i^i < oo) = nQl%s)A^i) = k\Ql%s)A^) = o, qI%s)A^i) > o)- 

(5) 

and finally we use this result to show that there is a p := p(A, f3{s),5) such that for G N, 

nZx,p{s),5{Di) = k\Ti = 0, < oo) = p(l - p)^ (6) 

3.1 The relationship between the epidemic process and the branching process 

Ball and Donnelly [7] proved that the epidemic process /^"^ (t) for n G N and the branching process 
Z(t) can be coupled in such a way that there exists a constant c > such that 

sup ^0 a.s. for n ^ oo. 

0<t<clog(n) 

Furthermore, observe that if d["'^ < oo, then /("^(t) > 1 for < t < and 
This, in turn, implies that for c > 0, 



["^ < t\D'^^ < oo) > 1 - e-'K 



lim P(Z)j"^ < clog(n)|L>f^ < oo) = 1. (7) 

n— >oo 

Combining the result by Ball and Donnelly with ^ gives that for any k €N, 
lim P(/(")(d|"^) = k\D\"^ < oo) 

n— >oo 

= lim P(/(")(L»i")) = < clog(n))P(L»i"^ < clog(n)|L'S"^ < oo) 

n — ^■oo 

+ lim P(/(")(Di"^) = k\D\''^ > clog(n))P(Di"^ > clogfn)!!?!"^ < oo) 

n^oo 

= lim P(/(")(d|"^) = k\D\"^ < clog(n)) 

n^oo 

= lim F{Z{Di) = fclZ^i"^ < clogn) 

n—yoo 

= ¥{Z{Di) = k\Di < oo). 

The above results allow us to analyse the branching process Z{t), instead of the SIR epidemic in 
large populations. So, from now on we will consider Z{t) instead of (t). 
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3.2 A random time change: from branching processes to the M/G/l-PS queue 

In this subsection we use a random time change to show that equation ([5|) holds. Let 

rt 



-{t) = [ l/Z{t')dt'. 
Jo 



We see that for « G N and T(t) < Ti, y4j(r(t)) = 0. Furthermore, 

^MT{t)) =1(0 < A{T{t)) < OO)^ =1(0 < A{T{t)) < ^)/Z{t), 

at at 

where Ti = and for i < j, Ti < Tj a.s. The random time change does not change the fact that the 
random variables Lj are i.i.d., distributed as L and independent of Tj for 1 < j < i. Furthermore, 

oo 

z{T{t)) = ^(i(r, < T{t)) -i{A{T{t)) = oo)), (8) 



i=l 



and 



P(r, G (r(t),T(t + h))\Tr^t^,Z{T{t)) = k^E^Tj < r{t)) = i - I) 

i=i 

= Xk{T{t + h) - T{t)) + o{h) (9) 
= Xk{h/k + o{h)) + o{h) 
= Xhl{k>0) + o{h). 

Similarly, we deduce 

oo 

P(A e (r(t),r(t + /i))|J^^(i),Z(t) = k, ^l{Dj <r(t)) = i - I) = 6hl{k >0)+o(/i), 

■'^^ (W) 
oo oo y^'^j 

P(El(7i e (r(t),r(t + /i))) + X:i(A e (r(t),r(t + /i))) > 1) = o{h). 

1=1 i=l 

Note that, as long as Z(r(i)) > 0, the description of the process -^a,/3(s),(5(''"(*)) is the same as 
the description of Qxp[s) 5(^)' '^ith Tj (resp. Di) replaced by Tj (resp. Di). So, 

nZx,(3is)ADi) = m = 0,^a,/3(.),5(^i) > 0) = P(qPS = ^f^ = 0, min g5^f(,),5(t) > 0). 

(11) 

We proceed by showing that 

P(Q(L>i) = A:|ri = 0, min Q(t) > 0) = P(Q(Di) = fc|Q(0) = 0, Q(Di) > 0). 

0<t<Z)i 

We use the following notation for the starting and stopping times of the busy periods of the queue: 

cri := min{t > 0;(3(t) > 0}, 

(7* := min{i > cj„; Q{t) = 0}, for n G N. 

Gn := m.m.{t > a*^^^; Q{t) > 0} forn E N \ {!}. 
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Furthermore, let N{t) := max{n > 1;(T* < t}, be the number of times the queue becomes empty 
in the interval (0, t). Observe that 

oo 

¥{Q{t + s) = k\ = s) = l,J^s) = P(Q(t) = k\Ti = 0), 

i=l 

and that Di is exponentially distributed and independent of the queue length upto time Di. There- 
fore, 

¥{Q{Di) = k\Q{0) = 0,Q(£»i) > 0) 
= E^=o nQiDi) = k\N{Di) = n, Q(0) = 0, Q{Di) > 0) 

P(iVpi) = n\Q{0) = 0, Q{Di) > 0) 
= En=onQiDi) = k\N{bi) = n,0 < Gn+i < Di, Q{Di) > 0) 

r{N{Di) = n\Q{0) = 0, Q{Di) > 0) 
= E7=oHQi[Di-'Tn+i]+^n+i) = k\N{Di) = n,0 = fi<an+i<Di,Q{Di)>0) (12) 

r{N{Di) = n\Q{0) = 0, Q{Di) > 0) 
= En=onQ{Di) = k\N{D,) = 0,ri = 0, Q{D,) > 0) 

P(iV(L'i) = n\Q{0) = 0, Q{Di) > 0) 
= F{Q{Di) = k\fi=0,Di<al) 
= ¥{Q{Di) = k\fi = 0,mino<,<^^ Q{t) > 0). 

Combining this result with equation (fTTjl leads to equation ([5]). 

3.3 The distribution of Q{Di) in the M/G/l-PS queue 
By [m eq. (2.6)], we know that for Q{t) = Q^^(^) 5, 

1 



e-%(s<3W|Q(o) = 0)dt 



5 + (l-s)A(l-7r) 



where vr is the smallest root of the equation vr = (3(6 + A(l — vr)) := E(e (^+-^(1 '^))^). Since 
^((x) := E(e~*^'^+''*(^~^))^) — X, is convex and g{0) > and g{l) < 0, vr is the unique root of g{x) = 
in [0, 1]. 

We observe that 

E(s^(^i)|Q(0) = 0) = 5 / e"'^*E(s^W|Q(0) = 0)dt 



5 + {l- s)A(l -^) 

and that if the random variable X is geometrically distributed with parameter p, then E(s^~^) = 
Yl'h=oP{^ ~ p)^s^ = p(l ~ (1 ~p)s)~^- By combining these observations with the fact that the 
probability generating function determines a distribution on the positive integers completely |15j . 
we deduce that 

P(Q(5i) = k\Q{0) = 0, Q(5i) > 0) = p{l - pf-^ (13) 

with 

p = 5/(5 + (i_vr)A), 

vr = {x G [0,l];a; = E(e-('^+(i-^)^)^))}, ^ ^ 
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i.e., conditioned on Q{0) =0 and Q{Di) > 0, Q{Di) is geometrically distributed (this has already 
been observed in [8]). 

Combined with the results of the previous subsection, we come to the main result of this paper: 
For ken 

IP(^A,/3(.),5pi) = k\D^ < oo) = p(l - pf-\ (15) 

with p as above. 

4 A special case: Markovian models 

The proof of the result by Kitaev [IT] discussed in l3.3l is rigorous, but we did not succeed in finding 
a intuitive argument for why it should be true. However, if the queueing process has the Markov 
property, then we can provide an intuitive proof. This will be done in this section. 

If L is exponentially distributed with parameter /i, then the SIR epidemic, the branching 
process and the queueing process under consideration are all Markovian. The branching process 
then becomes a simple birth-death process [I5l p. 251] with detections, described by the following 
equations: 

F{Z{t + h) = k + l\Z{t) =k) = Xkh + o{h), 
F{Z{t + h) = k-l\Z{t) =k) = fikh + o{h), 
F{Di&{t,t + h)\D,^i<t < Di,Z{t) =k) = 6kh + o{h), ^ ^ 

P(more than 1 event in {t,t + h)\Z{t) = k) = o{h). 

The M(A)/G(y^)/l-queue with exponentially distributed workloads with parameter //, is usu- 
ally referred to as the M{\)/M{fi)/1 queue [T5| p. 420], [9]. A property of this model is that many 
service disciplines lead to the same process Q(t), because no matter which customer is served, the 
served customer will leave the queue with probability fih + o{h) during an interval of length h. 
In particular, the First-In-First-Out (FIFO) discipline (in which customers are served 1 by 1 until 
they are fully served in the order of arrival), the Last-In-First-Out (LIFO) discipline (in which the 
server always serves the customer that arrived last out of the customers still in the queue) and the 
Processor Sharing discipline for M/M/1 queues with catastrophes all lead to the equations 

F{Q{t + h) = k + l\Qit) = k) = Xh + o{h), 

F{Q{t + h) = k - l\Q{t) = k) = nhl{k>0) + o{h), 

F{Di£{t,t + h)\Di^i<t<Di) = dh + oih), ^ ' 

P(more than 1 event in {t,t + h)) = o{h). 

We continue by considering an M/G/1 queue with catastrophes if the server uses a LIFO 
discipline, Q^^^^{t) = QY^^^^it). We deduce that for /c G N and some < p' < 1, 

P(qL^^O(^i) = k\bi < oo) = p'{l-p'f. 

It will turn out that p' = p, [8\. If L is exponentially distributed, then this result applies also to 
the processor sharing discipline, because as stated above, Q^^^^it) and (3^^(t) have the same law 
for M/M/1 queues. 
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4.1 The distribution of Q{Di) in the M/G/l-LIFO queue with catastrophes 

For the M/G/l-LIFO queue with catastrophes, the event {Ti = 0, mm Q{t) > 0} is the same 

0<t<Di 

as the event {Ti = 0, customer 1 is stih in the queue at time Di}. Furthermore, note that if a 
LIFO service discipHne is appUed and customer i is not in the queue at time Di, then none of the 
customers that arrived between the arrival and departure of customer i, wih be in the queue at time 
Di. So, we may ignore the arrival of such customers. Therefore, we say that the only two "good" 
events are arrivals of new customers, that will be in the queue at time Di and catastrophes. The 
first "good" event after Ti occurs during the period that customer 1 is served and the probabilities 
are (5/((5+ (1 — 7r')A) and (1 — 7r')A/((5 + (l — 7r')A) for a arrival and a catastrophe respectively, where 

tt' := P( min Q{t) = 0\ti = 0) 

0<t<Di 

is the probability that the queue is empty before the first catastrophe occurs. If this "good" event 
is a catastrophe, then Q[Di) = 1. If the first good event is an arrival, then customer 1 will not be 
served any more before -Di, and because Di is exponentially distributed, the number of individuals 
that arrived strictly after Ti and will add to Q{Di), is distributed as Q{Di). So, 

P(Q(£»i) = llTi = 0, min Q{t) > 0) = — — 

0<t<Di + (1 - TT'jA 

and for A; G N, 

P(g(L'i) = k + l\fi = 0, min Q{t) > 0) = ^ ^ , P(Q(-Di) = k\fi = 0, min Q{t) > 0). 

0<i<Di + (i - TT'jA 0<t<Di 

Combining this with the result of equation (|12p in section [3.21 gives that for k G N, ¥{Q{Di) = 
k\Q{0) = 0, Q{Di) > 0) = p'{l - p'f~^, where 

p' = 5/{5+{1-tt')\). (18) 

The same argument can be used to show that ¥{Q{Di) = 0|(5(0) = 0) = p' . Note that vr' is 
the probability that during the period customer 1 is served, neither catastrophes nor arrivals of 
customers that are still in the queue at Di, occur. Catastrophes occur at rate 5 and arrivals 
of customers that are still in the queue after an exponentially [5) distributed time occur at rate 
(1 — 7r')A. So, vr' is the smallest (and unique) root in [0, 1] of 

vr' = E(e-(^+(i--')^)^). (19) 

This is the same equation as the equation for vr given in the previous section and thus, 

¥{Q^^{Di) = fc|QPS(o) = 0) = P(QLIFOp^) ^ ^jgLIFO^Q) ^ o)_ 

In particular, we have provided a straightforward proof for the fact that P(Q(L)i) = A;|Q(0) = 
Q,Q{Di) > 0) = p{l - p)'', in the M/M/1 queue with catastrophes and ¥{Z{Di) = k\Di < oo) = 
p{l — p)^ in the simple birth-death process. 
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5 Possible applications 



5.1 Within herd spread of contagious animal diseases 

The spread of very contagious animal diseases like Classical Swine Fever (CSF), Foot and Mouth 
Disease (FMD) and Avian Influenza (AI) within farms usually progresses as follows: The pathogen 
is imported in a farm and one animal in the farm becomes infected. Then an SIR epidemic starts 
to spread within the farm. At a certain moment the farmer observes that an animal is diseased. 
Upon this detection the whole herd will be culled to prevent spread to other farms. In general the 
number of animals within a farm is large and the number of infectious animals at the moment of 
detection is small enough to accept the branching process approximation. 

In p] the following model for the spread of CSF within a farm is used. The spread starts 
with one infectious individual at time t = 0. Every infected animal "brings forth" new infected 
individuals at rate A. Infectious animals recover at rate n and are detected at rate 6. Upon 
detection all animals in the farm are culled. Because the animals are immediately culled upon 
detection, the distribution of the number of infectious individuals at the moment of detection is the 
best information we can hope for to obtain if we only consider the within-farm spread. However, 
questions on this distribution are not addressed in p]. 

After some straightforward algebra (114p leads to 

/i + (5 + A - J in + 5 + \y - 4A^ 
" = 2A 

and 

^ H + S + X- J{fi + (5 + A)2 - 4Xfi 

p = 1 . 

2/i 

Because the distribution of the number of infectious animals at the moment of detection is 
described by 1 parameter, more information than the number of infectious animals in a detected 
herd is needed in order to estimate r := X — n and R := X/n (respectively the exponential growth 
rate of the expected number of infected individuals and the basic reproduction number [S] within 
a farm). Possible further information should be obtained by contact tracing (by which we may 
obtain estimates for the real time between the first infection within a farm and the moment of 
detection) , or by looking for traces of immune response in all culled animals on a farm and obtain 
the total number of animals infected before detection. Note that this total number is geometrically 
distributed with parameter 5/{X + 6), because every time a detection, infection or recovery occurs, 
the probability that a detection occurs is 6/X times as high as the probability that an infection 
occurs. 

5.2 Spread of nosocomial pathogens like MRS A 

The spread of infectious diseases in hospitals is usually different from the spread of diseases in the 
general community. Because of antibiotic pressure and because of the weak immune responses of 
hospitalized patients, these people are more susceptible to many pathogens, than healthy people. 
In [2j the following model for the spread of MRSA in hospitals in countries with low prevalence of 
MRSA (like The Nordic countries and The Netherlands) is used: At very low rate MRSA carriers 
enter a hospital. Say that the MRSA is brought in at time t = and carriers arrive as single 
patients, so /(O) = 1. Because the rate at which MRSA carriers enter the hospital is very low. 
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we exclude further introductions from outside the hospital during the outbreak started by the first 
entering of an MRSA carrier. Every infected patient infects other patients at rate A, is discharged 
from the hospital at rate ^ and is detected at rate 6. As long as I{t) is small compared to the total 
number of patients in the hospital or ward, the constant infection rate is reasonable. 

Upon detection all patients in the hospital (or ward) are screened for MRSA colonisation and 
all infectious patients still in the hospital will be detected and isolated. Ideally this would make 
further infections within the hospital impossible. So I{t) increases by 1 at rate \I{t), decreases by 
1 at rate nl{t) and detection occurs at rate 51 {t). The dynamics of I{t) in this model are exactly 
the same as the dynamics of I{t) in the model for the within farm spread of contagious animal 
diseases. Therefore, all results of the previous subsection 15.11 can be used for the spread of low 
prevalence nosocomial pathogens. Note that the assumption that patients leave the hospital at a 
constant rate is not necessary for applying the results of this paper. 



5.3 Change of behaviour because of knowledge of the epidemic 

If an infectious disease is known to spread in a certain region, physicians will be more alert on 
symptoms of the disease and people will avoid crowded places or try to prevent their selves to 
become infected in other ways, like wearing masks as people did during the SARS epidemic in Asia. 
Therefore, it is reasonable to assume that the first detection of an infected person will lead to an 
increased rate of detection and a decreased rate of infection. 

Assume that the infection rate before the first detection was Ai/(t) and that the detection 
rate was 5il{t), while after the first detection the infection (resp. detection) rate will be X2l{t) 
(resp. 62! (t)). We assume that the recovery rate per individual, n does not change because 
of the first detection. From previous subsections we know that the number of infectious indi- 
viduals at the moment of the first detection is geometrically distributed with parameter pi := 

_ At+^i+Ai— y/(At+5i+Ai)2-4Ai/x 

Prom [15i p. 252] we deduce that the distribution of the number of individuals in a birth-death 
process at time r, /(r), with per capita birth rate A2 and death rate /i, which started with 1 
infectious individual is given by F{Z{t) = k) = qk{T) for k £ No, with 



ff2r _ I 

(1 - go(r))(l - R2qo{r)){R2qo{r)y-^ for all i > 1, 



where r2 = A2 — M and R2 = 

Some algebra yields that if the initial number of infectious individuals is geometrically dis- 
tributed with parameter pi, then at time r the number of infectious individuals Y{t) will be given 
by: 

p(y(r) = o) - 



1 - {l-pi)qo{T) 

P(y(r)=z) = (1 , J^'"'^^'°^:^^f\(l- (1 - forz>l 

^ ^ ' ^ l-{l-p,)qQ{r)'l-{l-pi)qo{Ty I - (I - p,)q,(r)' 

In this example an exponentially distributed infectious period is assumed, because then the 
number of infectious individuals at the time of first detection contains as much information for the 
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description of the progress of the epidemic after this detection, as knowledge of the whole process 
up to the time of first detection does. Note that for the epidemic with general infectious period it is 
possible to obtain the distribution of the number of infectious individuals that has been infectious 
for at least x time units at the moment of first detection from eq.(2.4)]. 

6 Extensions, limitations and questions 

In the previous sections we have assumed that the infection rate and detection rate during an 
infectious period are constant. It is tempting to conjecture that if the detection rate and infection 
rate are changing over time, but stay proportional, the number of infectious individuals at time Di 
is still geometrically distributed. Or formulated in the terminology of the corresponding branching 
process: let A(a) be the (possibly random) rate at which an individual at age a gives birth and 5{a) 
the (possibly random) detection rate of this individual at age a (if it is still alive at age a), where 
A(a) = c6{a) for some non-random constant c. 

An example of such a model is the SEIR (Susceptible Exposed — > Infectious Removed) 
epidemic model. In this model an individual first goes through a latent/exposed state after being 
infected, and after some random time, the individual becomes infectious itself. During the infectious 
period an individual can be detected, which happens at rate 5. Apart from the latent period the 
model is the same as the SIR epidemic model. 

If the infectious period is exponentially distributed with parameter and the latent periods 
are i.i.d. and distributed as A, then in the large population limit the dynamics of E{t) + I{t) are 
described by: 

¥{E{t + h) + I{t + h) = k + l\E{t) + I{t) = k, I{t) = I) = 
P(S(t + h) + I{t + h) = k- l\E{t) + I{t) = k, I{t) =1) = 
F{D^e{t,t + h)\Di^i<t<Di,I{t) = l) = 
P(more than 1 event in (t, t + h)\I{t) = I) = 

where events are infections, detections and recoveries. The time-change argument of Section [3] can 
be applied with 

r'(t) = / l{I{t') > 0)/I{t')dt'. 
Jo 

and we see that the dynamics of I{T'{t)) + E{T'{t)) are exactly the dynamics of an M/M/\ queue. 

However, simulations suggest that for the general SEIR epidemic, with non-exponentially dis- 
tributed infectious periods, E{Di) + I{Di) is not geometrically distributed in the large population 
limit. 

An open question is whether there are SEIR epidemics with non-exponentially distributed 
infectious periods, where E{Di) + I{Di) is geometrically distributed in the large population limit? 
Or if A(a) = c5{a) is a deterministic function, for which functions A(a) and which distributions of 
the infectious period L, I{Di) is geometrically distributed? To answer these questions it might be 
helpful to have intuitive understanding of why 

nQ{Di) = k\Q{Q) = 0, Q{bi) > 0) = p{l - pf-^ (21) 

holds for the general M/G/l-PS queue with catastrophes. The search for this intuitive understand- 
ing is still ongoing. 



\lh + o{h), 
filh + o(/i), 
6lh + o{h), 
o(h), 



(20) 
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